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Recently, there have been a number of fascinating and detailed experiments [|l], |2|, [| 
on layered antiferromagnets (AFMs) close to a zero-temperature (T) phase transition at 
which magnetic long-range-order (LRO) vanishes. The most prominent among these are the 
cuprates [|], |2| , which, upon doping with a small concentration of holes, lose their long-range 
Neel order and undergo a transition to a T = phase with magnetic, long-range spin-glass 
order; at a larger doping there is presumably a second transition to a quantum-disordered 
(QD) ground state. There have also been low-T experiments on layered AFMs on frustrated 
lattices ||, which have at most a small ordering moment. A remarkable feature of the 
measured dynamic susceptibilities of these AFMs is that overall frequency scale of the spin 
excitation spectrum is given simply by the absolute temperature. In particular, it appears to 
be independent of all microscopic energy scales e.g. an antiferromagnetic exchange constant. 

In this paper, we show that this anomalous dynamics is a very general property of finite 
T, 'quantum-critical' (QC) |4|] spin fluctuations near the initial onset of a T = QD phase. 
We present the first calculation of universal, QC dynamic scaling functions in 2+1 dimen- 
sions; these will be calculated for a model system - non-random, frustrated two-dimensional 
Heisenberg AFMs with a vector order-parameter. Quenched randomness will be shown to 
be a relevant perturbation to the clean system, and must be included in any comparison 
with experiments. Scaling forms for the dynamic susceptibility in random AFMs will be 
presented, and exponent (in-) equalities will be discussed. 

QC dynamic scaling functions can also be studied in other dimensions. In 1+1 dimensions 
the exact scaling functions can be obtained by a simple argument based on conformal invari- 
ance ||; most 3+1 dimensional models are in the upper- critical dimension and we expect the 
scaling functions to be free-field with logarithmic corrections. This leaves 2+1 dimensions, 
which is studied here for the first time in the context of AFMs: however, our results are more 



2 



general, and should also be applicable to other phenomena like the superconductor-insulator 
transition 

Most of our discussion will be in the context of the following quantum AFMs: 

H = Y^ JijSi ■ s i (!) 

where Sj are quantum spin-operators on the sites, i, of a two-dimensional lattice, and the 
Jij are a set of possibly random, short-range antiferromagnetic exchange interactions. The 
lightly-doped cuprates are insulating at T = 0, suggesting a model with completely localized 
holes: a specific form of TC with frustrating interactions was used by Gooding and Mailhot 
and yielded reasonable results on the doping dependence of the T = correlation length. 
Models with mobile holes have also been considered || and the results will be noted later. 

Two different classes of ground states of 7i can be distinguished: (i) states with magnetic 
LRO (Sj) = rrij and (ii) QD states which preserve spin-rotation invariance (Si) = 0. Further, 
we will distinguish between two different types of magnetic LRO: (A) Neel LRO in which 
case rrii ~ e 1 ® ' Ri with Q the Neel ordering wavevector and (B) spin-glass LRO in which 
case m, can have an arbitrary dependence on i, specific to the particular realization of the 
randomness. The lower critical dimension of the Heisenberg spin-glass Q may be larger than 
3 - in this case the spin-glass LRO will not survive to any finite T, even in the presence of 
a coupling between the layers. This, however, does not preclude the existence of spin-glass 
LRO at T = 0. 

Consider now a T = phase transition between the magnetic LRO and the QD phases, 
induced by varying a coupling constant g (dependent on the ratio's of the Jy in 7i) through 
a critical value g = g c , where there is a diverging correlation length £ ~ \g — g c \~ v ■ For Neel 
LRO, l/£ is the width of the peak in the spin structure-factor at the ordering wavevector 
Q. For spin-glass LRO, there is no narrowing of the structure-factor, and £ is instead a 
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correlation length associated with certain four-spin correlation functions JJ. At finite T we 
can define a thermal length £t ~ T~ x l z [z is the dynamic critical exponent) which is the 
scale at which deviations from T = behavior are first felt. The QC region is defined by the 
inequality £t < £ (Fig. 1); in this case the spin-system notices the finite value of T before 
becoming sensitive to the deviation of g from g c , and the dynamic spin correlations will be 
found to be remarkably universal. 

We consider first case (A) - a phase transition from Neel LRO to a QD phase; such transi- 
tions can occur for both random and non-random 7i. At T = the static spin susceptibility, 
X, will have a divergence at g — g c and wavevector q = Q: x(q — Qi 00 — 0) ~ \g ~ 9c\ ~ 7 
with 7 = (2 — 77)1/. The form of the q, u, T dependent susceptibility in the QC region can be 
obtained simply by finite-size scaling: £t acts a finite-size in the imaginary-time direction 
for the quantum system at its critical point j|, [| and hence implies the scaling form 

X(Q, w) = j^f® y Tl/Z , ^ J (2) 

where a\,d2 are non- universal constants, and $ is a universal, complex function of both 
arguments. The deviations from quantum-criticality lead to an additional dependence of $ 
on £t/£ : this number is small in the QC region and has been set to 0. Also of experimental 
interest is the local dynamic susceptibility Xl{w) — I dq x(<7>k>) = Xl + ^x'L with real 
(imaginary) part x'l (Xl) (Note dq = d 2 q/(4n 2 )). As x' ~ \q-Q\~ 2+v for \q-Q\ > cu^^T 1 ^, 
the real-part of the q integral is dominated by its singular piece only if 77 < 0. However x'l 
will involve only on-shell excitations, and the imaginary part of the q integral is expected to 
be convergent in the ultraviolet for both signs of 77. Thus the leading part of x'l wm always 
obey the scaling form 

with \i = rj/z, (4) 
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F(y) = y~^ / dx Im$(x, y) a universal function, and 03 a non-universal number, x'l a ^ so 
has a part obeying an identical scaling form which is dominant only if rj < 0. As we expect 
x'[ ~ w for small u;, we have the limiting forms F ~ sgn(y)|?/| 1_M for y <C 1 and F ~ sgn(y) 
for y> 1. Note that all the non- universal energy scales only appear in the prefactor 03 and 
the frequency scale in F is determined solely by T. 

Now the other case (B): the transition from spin-glass LRO to a QD phase. We do 
not expect singular behavior as a function of q because the spin-condensate rtij is a random 
function of i\ therefore the scaling form (Q) will not be obeyed. However the local suscep- 
tibility XlK) = Jo /{kBT) dre luJnT C(T), C(r) = (Si(0) ■ ftfrj) (where u n (r) is a Matsubara 
frequency (time) and the bar represents average over sites 1) will be quite sensitive to spin- 
glass LRO. In the spin-glass phase at T = lim r ^ 00 C(r) = mf > 0. In the QD phase, 
numerical studies of random, spin-1/2 AFMs [plj suggest that at T = 0, C(r) ~ l/r 1_Q , 



a > 0, for large r. At the critical point g = g c and T = we therefore expect the interme- 
diate scaling behavior with C(r) ~ l/r 1+M , x'l ~ Xl ~ an d — 1 < ^ < 0. In the QC 
region the scaling form @ for \l continues to be valid, despite the inapplicability of @. The 
limiting forms for F are as in (A), although the value \x is different: the Edwards- Anderson 
order-parameter obeys mf ~ \g — 9c\^ for g < g c - connecting the form of \l in the spin-glass 
phase to the critical point, we get 

H = -l + P/(zu) (5) 

We now consider various model systems for which exponents and/or scaling functions 
have be computed. 

We consider first a transition from Neel LRO to a QD phase in a non-random spin-1/2 



square lattice AFM with e.g. first (J x ) and second (J 2) neighbor interactions JIT, [12|]. As 



has been discussed in great detail elsewhere [12|, spin-Peierls order appears in the QD phase 



in this case (and in all other non-random AFMs with commensurate, collinear, Neel LRO). 
We now argue that the two-spin, QC dynamic scaling functions are not sensitive to the 
spin-Peierls fluctuations, and one may use an effective action for only the Neel order. It 
was found in the large M calculations for SU(M) AFMs that the asymptotic decay of the 
spin-Peierls correlations are governed by a scale £sp, which is much larger (for M large) 
than the scale, £, governing the decay of the Neel order The two scales are related 

by £ SP = £ 4M p\ where pi is a critical exponent given by p\ = 0.062296 to leading order in 
1/M D.S. Fisher |14| has noted that this is reminiscent of three-dimensional statistical 



models with a 'dangerously irrelevant' perturbation: e.g. the d = 3 classical XY model with 



a cubic anisotropy [15| has a phase-transition in the pure XY-class, but the 'irrelevant' cubic 
anisotropy becomes important in the low-T phase at distances larger than (■?/>> 1). By 
analogy, we may conclude that the spin-Peierls fluctuations are irrelevant at the critical fixed 
point governing the quantum phase transition, and relevant only at the strong-couping fixed 
point which governs the nature of the QD phase. 

It has been argued in Ref. that the dynamics of the Neel order-parameter is well 
described by an 0(3) non- linear sigma (NLa) model in the renormalized classical region 
(Fig. 1). The gist of the above arguments is that this mapping continues to be valid in the 
QC region - but not any further into the QD phase! We have computed properties of the 
QC region by a 1/N expansion on a O(N) NLa model: 

1 / dn\ ' 



Sn = — / drd 2 r 
2gJ 



(Vn) 2 + 



dr 



n 2 = 1 (6) 



were h is a real iV-component spin field, and c is a spin-wave velocity. The saddle-point 



equations of the large N expansion |16j were solved and the correlation functions were 
shown to satisfy the scaling forms (Q^J) to order 1/N (two- loops). We determined the 
values of $(x,?/) for real frequencies y by analytically continuing the Feynman graphs and 
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subsequently numerically evaluating the integrals. The numerical computations required the 
equivalent of 40 hours of vectorized supercomputer time. 

Our results for Im$ and F for N = 3 are summarized in Figs. 2,3. The transition has 
the exponent z = 1 which fixes the constant a-i = hc/ks in Eqn. (|2]). We normalized y) 
such that 9$ _1 /9x 2 |o,o — 1- Analytic forms for $ can be obtained in various regimes. We 
have 

Re$ _1 = Cq 2 + x 2 + ... x,y small (7) 
The universal number Cq, to order 1/N, is: 

Cq 1 = 0(1 + 0.22/AO ; 6 = 2 log ((1 + VE)/2 

Im$ has a singular behavior for x, y small: Im$(x = 0,y) ~ exp(— 3Q 2 / (2\y\)) / N while 
Im$(?/ < x) ~ yexp(— 30 2 /(2|x|))/iV\ With either x, y large, $ has the form 

^ = D Q {x 2 -y 2 Y l+ ^ 2 + ... ■ Dq = 1 — 0.3426/iV. (9) 



The exponents r/ have the known [17] expansion fM — rj — 8/(3n N) — 512/(277r N ) > 



The scaling function for the local susceptibility, F(y), has the limiting forms 

^) = sgu( 2 /)^M 1 -" y<l ; F(y) = sgn(y)^ Sin ^J 2 /2) y»l (10) 

As 77 is small, F is almost linear at small y. At iV = 00, F = sgn(y)9(\y\ — 0)/4. 

We now study Ti with quenched randomness. The simplest model adds a small fluctuation 
in the J\ bonds of the J1 — J2 model above: i.e. J± — > Ji+5J± where 5J\ is random, with r.m.s. 
variance J\, ensuring that a Neel LRO to QD transition will continue to occur. However 



the transition will not be described by the 'pure' fixed point as ^pure = 0.705±0.005 |Tg and 



thus violates the bound u > 2/d = 1 required of phase transitions in random systems [[19 



At long wavelengths we expect the spin fluctuations to be described by the NLa model, Sn, 



(§) with random, space-dependent, but time- independent, couplings g,c. A soft-spin version 
of Sfi with random couplings in d = 4 — e — e T space dimensions and e T time dimensions has 



been examined in a double expansion in e, e T [2C]. The expansion is poorly-behaved, and for 



the case of interest here (N = 3, e = 1, e T = 1) the random fixed-point has the exponent 
estimates r\ = —0.17, z = 1.21, v = 0.64, \i = —0.15. Note (i) fi,T] are negative, unlike the 
pure fixed point and (m) v is smaller than 2/d suggesting large higher-order corrections. 

Consider next a 7i on the square lattice with only J\ couplings, but with a small concen- 
tration of static, spinless holes on the vertices; this model will display a Neel LRO to QD 
transition at a critical concentration of holes. In the coherent-state path-integral formulation 
of the pure model, each spin contributes a Berry phase which is almost completely canceled 



in the continuum limit between the contributions of the two sublattices ||21|| . The model 
with holes will have large regions with unequal numbers of spins on the two sublattices: 
such regions will contribute a Berry phase which will almost certainly be relevant at long 



distances. Therefore the field theory of Ref. [20] is not expected to describe the Neel LRO 



to QD transition in this case. A cluster expansion in the concentration of spins has recently 



been carried out by Wan et. al. f22| : and yields the exponents i] = —0.6, z = 1.7, v = 0.8, 
\i = —0.35. Note again that /x, r\ < 0, although the violation of v > 2/d suggests problems 
with the series extrapolations. 

Finally, we have also considered || the consequences of mobile holes in a non-random 



AFM. The spin-waves and holes were described by the Shraiman-Siggia |23j field theory. 
Integrating out the fermionic holes, led to a spin- wave self-energy ~ a>i\q— Q\ 2 + a2UJ^ + . . ., 
(ai,02 constants) at g = g c , T = 0; non-analytic \u n \ terms appear only with higher powers 
of \q — Q\,u n indicating that the Neel LRO to QD transition has the same leading critical 
behavior as that in the undoped, non-random J\ — J2 model above. The exponents and 
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scaling functions are identical, but the corrections to scaling are different. 

To conclude, we discuss implications for neutron scattering experiments in the doped 
cuprates 0. The significant low T region with a T-independent width of the spin structure- 
factor indicates that the experiments can only be in the QC region (Fig. 1) of a T = 
transition from spin-glass LRO to QD: the diverging spin-glass correlation length will then 
not be apparent in the two-spin correlations. The numerical results of Ref. |7j also indicate 
that, in the absence of a coupling between the planes, a spin-glass phase will appear at 
any non-zero doping. The experimental x'l nas been fit to a form I(\u\)F(%uj /ksT) |l[ 
which is compatible with the theoretical QC result (Q) if / ~ M^- A fit to this form for 
/ in Lai^Sro^CuO^ yielded /i = —0.41 ± 0.05 with all the predited points within the 



experimental error-bars [24]. As it appears that only random models have \i < 0, it is 
clear that the effects of randomness are experimentally crucial, confirming the theoretical 
prediction of their relevance. Further theoretical work on the QC dynamics of random 
quantum spin models is clearly called for. 

We thank B. Keimer, G. Aeppli, R.N. Bhatt, D.S., M.E. and M.P.A. Fisher, B.I. Halperin, 
D. Huse, N. Read, and A. P. Young for useful discussions. This research was supported by 
NSF Grants No. DMR 8857228 and PHY89-04035, and the A.P Sloan Foundation. 
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Figure Captions 

1. Phase diagram of 7i (after Ref. 0). The magnetic LRO can be either spin-glass or 
Neel, and is present only at T = 0. The boundaries of the QC region are T ~ \g — g c \ zv ■ 
For non-random Ji which have commensurate, collinear, Neel LRO for g < g c , all of 
the QD region (g > g c ) has spin-Peierls order at T = 0-this order extends to part of 
the QD region at finite T. 

2. The imaginary part of the universal susceptibility in the QC region, <3>, as a function 
of x = hcq/ \ksT) and y = Tioj/(kBT) for a non-random square lattice AFM which 
undergoes a T = transition from Neel LRO to a QD phase. The results have been 
computed in a 1/N expansion to order 1/N and evaluated for N = 3. The two- 
loop diagrams were analytically continued to real frequencies and the integrals then 
evaluated numerically. The shoulder on the peaks is due to a threshold towards three 
spin-wave decay. 

3. The imaginary part of the universal local susceptibility, F, for the same model as in 
Fig 2. We have F(y) = y^^ J dx Im$(x, y). The oscillations at large y are due to a 
finite step-size in the momentum integrations. 
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